Nodal Superconducting Order Parameter and Thermodynamic Phase Diagram of (TMTSF) 2 C10 4 
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The organic materials (TMTSF^X are unique unconventional superconductors with archetypal quasi-one- 
dimensional (Q1D) electronic structures. Here, based on our comprehensive field-angle-resolved calorimetry 
of (TMTSF)2C1C>4, we succeeded in mapping the nodal gap structure for the first time in Q1D systems, by 
discriminating between the Fermi wavevectors and Fermi velocities. In addition, the thermodynamic phase 
diagrams of (TMTSF)2C104 for all principal field directions are obtained. These findings, providing strong 
evidence of nodal spin-singlet superconductivity, serves as solid bases for further elucidation of anomalous 
superconducting phenomena in (TMTSF^X. 



It was in 1980 when the first organic superconductor 
tetramethyl-tetraselena-fulvalene (TMTSF) salt was discov- 
ered. [1, 2] Since then, (TMTSF) 2 X (X = C10 4 , PF 6 , etc.) 
have been widely studied because of their fascinating prop- 
erties resulting from their archetypal quasi-one-dimensional 
(Q1D) conductivity and strong electron-electron interac- 
tions. [3-6] Interestingly, the superconducting (SC) phase is 
located next to a magnetic phase, resembling superconductiv- 
ity in other unconventional superconductors such as high-T c 
cuprates and pnictides. [7, 8] Because of this similarity, as 
well as the simplicity of the electronic structure, investiga- 
tions of (TMTSF^X can provide useful guidelines for studies 
of other unconventional superconductors. 

Many studies have revealed unusual SC phenomena in 
(TMTSF^X. For example, the onset of superconductivity in 
resistivity is observed even above 4 T when the field is par- 
allel to the a axis, the most conducting direction, or to the b' 
axis, the second-most conducting direction. [9-12] This fact 
indicates that a certain contribution of superconductivity sur- 
vives beyond the Pauli limiting field fioHp ~ 2.3-2.6 T, [12] 
where ordinary singlet pairs would be unstable due to the Zee- 
man splitting. Thus, possibilities of a spin-triplet pairing state 
or a spatially-modulated spin-singlet pairing state, which is 
the so-called Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, 
have been discussed. [9-14] In addition, the high-field super- 
conductivity is accompanied by a peculiar anisotropy of the 
resistivity onset. [11, 12] However, SC phase diagrams have 
not been established from thermodynamic measurements. 

For clarification of the origin of such unusual behavior, as 
well as for identification of the SC mechanism of (TMTSF^X, 
orbital and spin parts of the SC order parameter is essentially 
important. Although sign changes on the SC gap is evidenced 
by the suppression of superconductivity by a tiny amount of 
non-magnetic impurities, [15, 16] details of the gap structure 
and the SC symmetry are still controversial. [17-19] Theories 
based on a simple Q1D model with the Fermi surface (FS) 
consisting of a pair of warped sheets have revealed that a spin- 
singlet (i-wave-like state with line nodes is stable when spin 
fluctuations drive the pairing, [4, 20] whereas a spin-triplet 
/-wave-like state with similar line nodes can be stable when 
charge fluctuations are incorporated. [4] The CIO4 salt have 



a slightly different FS due to the orientational order of the 
tetrahedral CIO4 anions at T\o = 24 K; this order leads to 
a folding of the band structure along the b* direction and the 
FS splits into two pairs. [21] A nodeless c/-wave-like, a node- 
less /-wave-like, and a nodal c/-wave-like states have been 
proposed for such a FS. [22, 23] 

To clarify these issues, we performed field-angle-resolved 
calorimetry for one piece of a (TMTSF^CICU single crys- 
tal. We used a crystal weighing as low as 76 \xg grown by 
an electrocrystallization technique. [24] This piece of crystal 
was previously used in our previous transport study, and was 
confirmed to be very clean with the mean free path of as large 
as 1.6 |im. [11, 12] We used a 3 He- 4 He dilution refrigerator to 
cool the sample. After we cooled the cryostat and sample to 
4.2 K, we heated the sample again to 26 K, and cooled it very 
slowly across Tao to 20 K at 4 mK/min, so that the anions 
order well to be in the "relaxed" state. We newly developed 
a high-resolution (~ 100 pJ/K at 1 K) calorimeter shown in 
Fig. 1(b), based on a modification of the "bath-modulating 
method". [25] The advantage of this technique is that a heater 
on the sample holder is not necessary; thus the background 
heat capacity of the sample holder can be minimized. Al- 
though this background contribution is not subtracted from 
the data shown here, we have checked that the background is 
nearly field-independent. [24] The magnetic field is applied 
using the vector magnet system. [26] The magnetic field is 
aligned to the crystalline axes by making use of the anisotropy 
in the upper critical field H C 2- The precision and accuracy of 
the field alignment is approximately 0.1°. More details of the 
experimental procedure will be described elsewhere. [27] 

First we focus on the field dependence of C/T plotted in 
Fig. 1. For H \\ c*, i.e. fields perpendicular to the conduct- 
ing ab' plane, C(H)/T exhibit a H 5 dependence at low tem- 
peratures. This dependence provides strong evidence for a 
line-node gap. [28] In contrast, C(H)/T for H \\ a, for which 
the orbital pair breaking is substantially weakened, exhibits a 
concave-up curvature at low temperatures (Fig. 1(a)) near H C 2- 
Considering the clear change of the spin susceptibility in the 
SC state, [18] we identify this behavior as the Pauli limiting 
behavior in a spin-singlet superconductor. [29] 

To investigate the nodal structure of the SC gap, we mea- 
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FIG. 1. (color online) Magnetic field dependence of C/T at 0.11 K 
for (a) H || a (red squares), H \\ V (blue circles) and (c) H \\ c*. The 
arrows indicate the onset H C 2- The broken curve in (b) illustrates H 05 
behavior. The inset in (b) is a photo of our calorimeter. 



sured the in-plane field-angle dependence of C/T represented 
in Fig. 2. When the SC gap has a node or a zero at fc p ode , 
the quasiparticle (QP) density of states (QDOS) N, which 
is proportional to C e /T at low temperatures, varies with the 
field direction. Here, C e is the electronic heat capacity. In 
the low-temperature and low-field limit, such variation of 
QDOS originates from the "Doppler shift" of the QP energy 
5u)(r, k) oc v s (r) ■ Vp(k), where v s (r) is the velocity of the 
supercurrent around a vortex and Vp(k) is the Fermi veloc- 
ity. [28, 30] When 5a> is larger than the gap A(fc), QPs with 
the wavevector k are excited. Because most of the excitation 
occurs in the vicinity of the nodes, the most important term 
is c5w node oc v s ■ t>™ de , where t>™ de = t> F (fc F ode ). When H is 
parallel to v™ de , 6u) node becomes zero (i.e. v s _L v™ de ) be- 
cause v s is perpendicular to H, and QDOS induced by this 
shift becomes small. (Fig. 4(a)) Using this idea, we can in- 
vestigate the directions of Vp at nodal positions. If the con- 
dition kp || vp, is satisfied, the £-space nodal direction equals 
the field direction for which C/T exhibits minimum. [31] For 
Q2D or 3D systems, the assumption kp \\v? is reasonable as a 
simple model. In contrast, for Q1D systems, Vp is not parallel 
to feF even in the simplest model. Thus, to deduce the nodal 
position in the £-space, information on the band structure is 
required. Another difficulty for Q1D systems is that a large 
in-plane anisotropy in H C 2 also contributes to the anisotropy in 
C e . Because of these difficulties, C e /T oscillation in Q1D sys- 
tems has been little studied, despite its essential importance in 
determining the SC gap structure. 

Interestingly, we find that the C{</>)/T curves of 
(TMTSF)2C104, where <p is the azimuthal angle of the field 
measured from the a axis, become asymmetric with respect 
to the a axis (e.g. C((f>) > C(-0) for 0° < < 90°) at low 
temperatures and low fields, as shown in Figs. 2(a) and (b). In 
contrast, the curves are nearly symmetric at high temperatures 
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FIG. 2. (color online) Magnetic-field-angle <p dependence of the heat 
capacity for fields rotated within the conducting ab' plane at (a)-(c) 
0.14 K, and (d)-(f) 0.50 K. For comparison, the same data are also 
plotted against -<p with appropriate shifting (the black curves). The 
deviation of the two curves indicates the asymmetry in the C((f>)/T 
curve, (g) QDOS iV(0)/iV(9O°) (blue solid curve) and N(-<f>)/N(90°) 
(black broken curve) calculated by eq. (1) assuming two nodes (nl 
and nl) with the parameters 0„i = -10° and 0„2 = +10°, A„2/A n i = 
0.3, and T = // c2 (0°)/// c2 (90°) = 3.5. 



or in high fields. Thus the asymmetry is not due to a misalign- 
ment of the field. What is more, at 0.14 K, small kinks are 
observed at around (p = +10° as shown in Fig. 3(a). The kink 
signatures are more obvious in the derivatives (Fig. 3(b)): The 
step-like behavior in the first derivative and the peaks in the 
second derivative provide clear indication of the kinks. 

The results qualitatively agrees with the theoretical expec- 
tation that the specific-heat anomaly due to the gap anisotropy 
should appear only in the low-temperature and low-field 
region [31, 32]. Thus we attribute them to the SC gap 
anisotropy. Below, we demonstrate that a simple model based 
on the Doppler shift reproduces the key features of the un- 
conventional behavior. In order to incorporate the large in- 
plane H C 2 anisotropy into the Doppler-shift mechanism, [30] 
we assume that QDOS N{<p) for H «c H C 2(</>) varies as 
N((f>) oc ^H/H C 2(</>) Yin An\ sin(0 - 0„)|, where <p„ is the di- 
rection of vp at the «-th node/zero. To take into account the 
triclinic band structure, we introduce the nodal-position de- 
pendent amplitude A n . We approximate // C 2(0) to follow the 
effective mass model: H c2 ((f>) = H c2 (0°)/(r 2 sin 2 0+cos 2 0) 1/2 , 
with T = H C 2(0°)/H c2 (90°). Thus N((f>), which is proportional 
to C e /T, should vary as 
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FIG. 3. (color online) (a) Enlarged view of C(<p)IT at 0.14 K and 
0.3 T. The arrows indicate the positions of the small kinks, (b) 
Derivatives d(C/T)/d0 (green crosses) and d 2 (C/T)/d(p 2 (red trian- 
gles), (c) Calculated DOS near <p = 0. (d) Qualitative behavior of 
dN/d(p (green broken curve) and d 2 N/d(f> 2 (red solid curve). 



The calculated N(<f>) with <f>„ = ±10° plotted in Figs. 2(g) 
and 3(c) well captures the observed unconventional behav- 
ior, though we could not perfectly fit the equation to the data. 
For a better fitting, additional contributions of vortex cores, 
thermal excitations, and the Pauli effect should probably be 
taken into account. It is worth noting that the behavior of the 
derivatives is also consistent with the experiment. We note 
that a recent theory [32] also attributed the observed calori- 
metric behavior to the existence of the gap node. 

The above analysis indicates that at least there must be 
a node (or zero) with v$ pointing (f> = +10° and another 
one with v F pointing -10°. To determine the nodal po- 
sition in the £-space, we plot in Fig. 4(b) the k y depen- 
dence of the velocity angle measured from the a axis, <p Vf , 
based on the tight-binding band structure [6, 21]. Fig. 4(b) 
manifests that \<j> v? \ reaches 10° on the outer FS at k y ~ 
±0.25fc*, +0.36/7*, -0.06b*, where b* is the size of the 1st Bril- 
louin zone along Ic y . Thus, at least, some of the gap nodes 
or zeros should be located around these positions. The sim- 
plest structure that satisfies this condition is a structure with 
nodes or zeros running at ky = +0.25b*, which is presented in 
Fig. 4(c). Leading candidates of the pairing state that satisfies 
the obtained gap structure are the cZ-wave-like or the g-wave 
like states shown in Fig. 4(c). 

Finally, in Fig. 5, we present the thermodynamic SC phase 



FIG. 4. (color online) (a) Quasiparticle excitation due to the Doppler 
shift. When the field is parallel to v F at a node, the quasiparticle exci- 
tation at this node is reduced, (b) Dependence of 0„ F = arctan(Vy/v Jt ) 
on k y for the outer FS (O-FS; solid curve) and inner FS (I-FS; broken 
curve). The arrows indicate the points where </>„ F reaches ±10°. (c) 
Most plausible gap structure with nodes or zeros at ky ~ ±0.25b*. Ex- 
amples of the superconducting states that satisfy the observed nodal 
structure are shown: the li-wave-like state and the ^-wave-like state. 



diagrams based on the calorimetry compared with transport 
phase diagrams based on the c*-axis resistivity. [11, 12] The 
thermodynamic upper critical field uoH'~ m (0) 2.5 T for H || 
a is much smaller than noH"™ = -0.73T c (fiodH C 2/dT\T=Tj ~ 
7.7 T expected for the orbital pair-breaking but agrees with 
the Pauli-limiting field fioHp ~ 2.3-2.6 T. [12] This fact again 
supports a spin-singlet scenario. Furthermore, the absence of 
multiple SC phases provides the first proof of a singlet state 
in the region below //* m . Looking again at earlier thermody- 
namic studies here, one notices that the present ff ^(O) agrees 
with the field at which the nuclear-lattice relaxation rate re- 
covers to the normal-state value, [18] or with the irreversible 
field in the torque measurement. [10] Thus, these anomalies 
reported earlier are now turned out to be due to the thermody- 
namic SC transition. 

Although we used the identical crystal in both the thermo- 
dynamic and transport studies, the obtained phase diagrams 
are quite different. It is clear that a long-range ordered SC 



state does not exist above H^ m - In the region between 
and the resistance onset, however, the observed sharp resis- 
tance drop and the anisotropy in the onset temperature reflect- 
ing the bulk Fermi surface anisotropy [11, 12] evidence that 
intrinsic superconductivity robustly contributes to the trans- 
port. The superconductivity above ffS™ should be realized 
without any noticeable entropy change compared to the nor- 
mal state. Thus, it must be characterized either by a fluctuat- 
ing order parameter or by a static order parameter accompa- 
nied by only a tiny change in the density of states. The true 
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FIG. 5. (color online) Thermodynamic SC phase diagrams of 
(TMTSF) 2 C104 compared with the transport phase diagrams for 
fields parallel to (a) the a axis, (b) the b' axis, and (c) the c* axis. The 
closed symbols are obtained from field sweeps and the open sym- 
bols from temperature sweeps. For comparison, we also present the 
onset temperature of the c*-axis resistance R c , (open diamonds), as 
well as the temperature where R c . becomes zero (crosses), reported 
in Ref. 11. The broken lines indicate the slope po&H € 2(T)IAT near 
T = T c : -8.1 T/K for H || a, -2.3 T/K for H \\ b', and -0. 1 1 T/K for 
H || c*. The crystalline structure and axes are schematically shown. 



nature of the superconductivity in the high-field region is one 
of the most intriguing problems to be clarified in future. 

In summary, our precise field-angle-resolved calorimetry 
using a very clean single crystal of (TMTSF^CICU provide 
strong evidence of a spin-singlet nodal superconductivity. Our 
result marks the first calorimetric mapping of the gap structure 
of a Q 1 D superconductor beyond the constraint previously be- 
lieved; this technique is indeed applicable regardless of crys- 
talline and electronic symmetries. We also revealed the ther- 
modynamic SC phase diagram, which exhibits substantial de- 
viation from the resistivity onset. The present findings should 
serve as solid bases for investigation of the pairing mechanism 
of (TMTSF^X, as well as for elucidation of unconventional 
SC phenomena in this compound. 
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Supplemental Material 



SAMPLE 



In the present study, we used a single crystal weighing as low as 76 u.g shown in Fig. 6(a), whose heat capacity at temperatures 
slightly above T c is estimated to be approximately 3 nJ/K based on the result in the previous study. [33] It should be noted that 
this piece of crystal was previously used in our previous transport study [11]; we removed the electrodes from the crystal and 
used it for the present calorimetry. This allows us to compare superconducting properties of the identical sample. A small 
amount of gold used for the electrode is left on the surface of the crystal as seen in the photo. The contribution of this gold to 
the heat capacity is estimated to be on the order of 1 pJ/K at 1 K and is negligible. 

We note here that (TMTSF^X salts are particularly stable in the air and at room temperature over many years and don't suffer 
from aging. Indeed, recently published transport work has been performed on excellent 30 years old crystals. [7] Thus, the 
oldness of our sample is not a problem, although a few years have passed since we used the crystal for the previous transport 
study. We also note that thermal cycling of (TMTSF^X crystals are not reported to affect bulk measurements, although transport 
measurements are known to be affected by temperature cycling mainly by micro-cracks within crystals. 



TEMPERATURE DEPENDENCE OF THE HEAT CAPACITY AND BACKGROUND 



Our high-resolution (~ 100 pJ/K at 1 K) calorimeter used in this study is shown in Fig. 6(a). This calorimeter is developed 
based on a modification of the "bath-modulating method": [25] For our apparatus, small thick-film resistors are used for ther- 
mometers, though thermocouples are used in Ref. [25]. The advantage of the bath-modulating technique is that one does not 
need to place a heater on the sample stage; thus the background heat capacity of the sample holder can be minimized. With this 
technique, sensitivity high enough to measure the heat capacity of one single crystal of (TMTSF^CIC^ is achieved, although 
we can only obtain the relative value of the heat capacity. Depending on the temperature and field conditions, we chose the 
temperature modulation frequency as 34 Hz or 10 Hz, which is twice of the heater-current frequency. We adopted the temper- 
ature modulation amplitude of 1.5%-rms and 3%-rms of the base temperature for the temperature and field sweeps, while we 
used 6%-rms for the field-angle sweeps in order to improve the signal-to-noise ratio. More details of the measurement will be 
described elsewhere. 

The heat-capacity data presented in this report contains the contributions from both the sample and the background. The 
background contribution, which is represented in Fig. 6(b), is dominated by the heat capacity of the thermometer and that of 
the Apiezon N Grease used to fix the sample to the calorimeter. Because our single crystal weighs only 76 [ig, the background 
contribution is affected even by a tiny amount of the grease. The amount of the grease used in the measurement cannot be exactly 
known. Therefore, the subtraction of the background contribution is not thoroughly accurate. 




FIG. 6. Descriptions of our experimental setup, (a) Photos of the calorimeter and the sample. We used a very clean single crystal of 
(TMTSF)2C104 weighing 76 u.g. We note that this crystal is identical to the one used in our previous transport studies. [11, 12] (b) Temperature 
dependence of the total (i.e. sample + background) heat capacity (solid circles), compared with the background heat capacity (broken curve). 
Note that the background contribution is not thoroughly accurate as described in the text. 



FIELD DEPENDENCE OF THE BACKGROUND 



In order to prove that the background contribution have little influence in the field-amplitude sweeps and field-angle sweeps, 
we compare the data with the separately measured background data in Fig. 7. Although the actual background may be different 
by a few tens of percent, its field dependence should be quite similar to the data presented here. As is clear in this figure, the 
field-amplitude and field-direction dependences of the background contribution are negligibly small. Therefore, it is certain that 
our conclusion of the Letter is not at all affected by the background contribution. 
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FIG. 7. Comparison of the experimental data with the separately measured background. In all cases, the background does not affect the field 
dependence of the total heat capacity at all. (a)-(b) Field-sweep data for H \\ a and H \\ b' at 0.90 K and 0.1 1 K. (c)-(d) Field-sweep data for 
H || c* at 0.90 K and 0.1 1 K. (e)-(f) In-plane field-direction <p dependence at 0.50 K and 0.14 K. 
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